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1. Introduction 

The general feature of asymptotic freedom makes weak coupling approaches very natural in 
non-abelian gauge theories, such as the SU(3) model, which describes the gluonic degrees of free- 
dom of Quantum Chromodynamics. At asymptotically high temperatures low orders of perturba- 
tion theory may be acceptable, but at any lower scale that could be probed by a realistic experiment 
an extension is necessary: either by the inclusion of very high order diagrams, or by an efficient 
resummation scheme, such as Hard Thermal Loops (HTL). Note however that analytic perturba- 
tive expansions are plagued by infrared divergences due to which the series can be computed only 
up to a given finite order. There is strong simulation evidence that at low temperatures {Tc ~ 260 
MeV) the gluonic matter freezes and a first order transition takes place. At even lower temper- 
atures colorless non-perturbative excitations govern the thermodynamics. To describe the phase 
transition or the glueball gas no weak coupling scheme succeeds and one has to rely on a natively 
non-perturbative approach, such as lattice field theory. 

The past year witnessed considerable achievements on the side of the analytical results. The 
HTL scheme has been used to calculate the pure SU(3) gauge theory's thermodynamic potential to 
the next-to-next-to-leading order (NNLO) The authors used their results at intermediate tem- 
peratures {r^ ATc) where existing lattice data were available. Later the same authors have extended 
their results to full QCD (with massless quarks) |^ and found good agreement with the lattice data 
of the Wuppertal-Budapest collaboration from about 300 MeV. 

In conventional perturbation theory even higher orders can be computed by dimensional re- 
duction. The full expression up to g^log{g) order is given in [Q] and was compared to the Bielefeld 
lattice data at T = 4.57^. Fitting the pressure (thermodynamic potential) the slope of the pres- 
sure curve was succesfuUy predicted. This raised hope that at this high order perturbation theory 
does possess some predictive power at phenomenological temperatures. In this work we repeat this 
fitting procedure at a much higher temperature, where the sixth order can be shown to be a minor 
correction to the fifth order. Instead of treating the soft sector strictly perturbatively a screened per- 
turbation theroy can be formulated for the dimensionally reduced theory, resulting in significantly 
better convergence of the free energy 

For more than a decade the renowned paper by Boyd et al [||] has been the reference lattice 
simulation of the SU(3) theory in the temperature range of 1 . . .4.57^. It uses the plaquette gauge 
action at up to Nt = ^ lattice spacing and an aspect ratio of 4. Here Nt denotes the number of lattice 
points in the Euclidean time direction, meaning that the lattice spacing at any given temperature 
r is a = \/{TNt). The fixed Nt approach has been introduced in Ref. and this work follows 
it, too. It implies that the lattice spacing varies with temperature. Continuum limit is achieved 
by performing an 1 /Nt — )• extrapolation on the data at a set of fixed physical temperatures. The 
aspect ratio r = LT sets the ratio between space and time-like lattice points. 

Since the publication of [^] several similar simulations were performed to study pure gauge 
theory. The equation of state has been recalculated using the Symanzik improved gauge action [||]. 
This set of simulations have been further generalized to SU(A'^,.) theories with Nc > 3 in Refs. p^ ]. 
Alternatively, the equation of state can also be calculated by fixing the lattice spacing, and using Nt 
for tuning the temperature [pT]|. This approach is mostly advantageous with Wilson-type dynamical 
fermions, and less economic for the pure gluonic theory. 



2 



Lattice SU(3) thermodynamics and the onset of perturbative behaviour 



Sz. Borsanyi 



In most fixed Nt simulation projects, like Ref. [Q], the aspect ratio is kept constant to allow 
the use of a single lattice geometry. This means that higher temperatures are simulated at smaller 
volumes. Length scales that are present in (resummed) perturbative calculations, such as ~ T , 
~ gT and ~ g^T are normally well accommodated in the lattice, since the renormalized coupling 
g drops only logarithmically as the temperature increases. Yet, to establish the range of validity of 
the perturbative approach itself, one has to simulate the non-perturbative ~ Tc scale, too. In this 
sense, the aspect ratio sets the maximum temperature as a precondition for the non-perturbativeness 
of the simulation: T < rTc. In most previous works this was set to r = 4. 



2. Simulation setup 

In this work we calculate the continuum equation of state of the SU(3) theory using tree-level 
Symanzik improvement in the temperature range of T /Tc = 0.7 ... 16 (on 80^ x 5, 96^ x 6 and 
1 14^ X 7 lattices). We also support this continuum extrapolation using an additional Nt = S set of 
lattices (64^ x 8) below STc. Furthermore we present a non-continuum data set that is valid up to 
approximately 24 (on a 120'' x 5 lattice) and study finite volume effects using various smaller 
boxes. Finally, from our third set of simulations we calculate the continuum equation of state in a 
small box (on 40^ x 5, 48^ x 6 and 64^ x 8 lattices) up to 1000 Tc. These latter data we use to find 
the optimal free parameters of existing perturbative calculations, e.g. the preferred renormalization 
scale for the HTL scheme. The precision of our data points exceeds any previous calculation by 
about an order of magnitude. 

Although the lattice techniques for the analysis as well as for the generation of lattice con- 
figurations have been well established, achieving the presented statistics was a very challenging 
procedure. As explained in Ref. [^, the normalized trace anomaly {e — 3p)/T'^ is determined 
first. The normalized pressure p/T^ and energy density e/T^ are calculated using thermodynamic 
relations. The trace anomaly contains a quartic divergence, which is subtracted using lattice mea- 
surements of the same quantity (at the same lattice spacing) at a smaller temperature. Because of 
this divergence the trace anomaly is rather difficult to measure. Moreover, the value of {e — 3p)/T^ 
reduces rapidly with temperature as one moves away from the transition region in either direction. 

Another challenging issue was the accurate determination of the non-perturbative beta func- 
tion corresponding to the Symanzik improved action. Instead of using the string tension or the 
Sommer parameter, it was advantageous to define the lattice spacing in terms of the transition 
temperature. To this end we determined the critical couplings up to Nt = 20 from the peak of the 
Polyakov loop susceptibility. Matching to the universal two-loop running (in terms of the improved 
coupling in the "E" scheme [|l^], generalized for the case of the Symanzik improved action) we 
determined the lambda parameter in terms of the transition temperature: Tc/A-^ = 1.26(7). (The 
error is overwhelmingly systematic and reflects the sensitivity to various continuum extrapola- 
tions.) This is consistent with the combination of previous determinations: the Lambda parameter 
^MS ~ 0.614(2)(5)rQ ' of [13] can be translated to ^/a units using ^/oro = 1.192(10) (based 
on [14]) and then used with Td \fo = 0.629(3) of Ref. [||]. Through our direct result for Tc/Aj^ 
one can easily translate the scale setting of the perturbative expression for the SU(3) free energy to 
the lattice language. 
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Given the beta function one can simply relate the expectation value of the gauge action Sg (a 
weighted sum of 1 x 1 and 1x2 plaquettes) to £ — 3/? [0]. In the standard lattice renormaliza- 
tion scheme this {Sg)T result is then subtracted from the value corresponding to zero temperature. 
This would require a lattice with a large temporal size. In order to be able to fit larger lattices 
into the given memory of our computer system we used half-temperature subtraction, i.e. we cal- 
culated {e{T) - 3p{T)) /T"^ - ^{e{T /2) - 'ip{T /2) / {T /if along the lines of our previous work 
[15]. To finally arrive at (e(r) — 3p{T))/T^ this partial result was supplemented with another set 
of simulations at half temperature, double lattice spacing, but same physical volume. For these 
supplementary data we used the standard renormalization subtracting the vacuum. The continuum 
limit from this combined technique is equal to what one finds using the stardard scheme. 



3. Results 

We start the presentation of the results with the reproduction of the Boyd et. al. data [^] in the 
left side of figure [!} We fit our four large- volume data sets with different lattice spacings altogether, 
using an A^,-dependent spline function. As a result we have a smooth function interpolating our 
data for each Nf (colored lines in the figure), together with a smooth, continuum extrapolated curve 
(yellow band in the figure). The systematic error coming from this extrapolation procedure and 
the statistical error are added in quadrature. We see a small discrepancy between our results and 
Ref. [Q], which might be attributed to the scale setting assuptions in [|5|]. 
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Figure 1: Left side. The trace anomaly calculated on 80^ x 5, 96^ x 6, 114-^ x 7 and 64^ x 8 lattices. 
The continuum extrapolation makes use of the N,''^ scaling of the Symanzik action. The data deviates from 
Boyd et al |Q especially close and below the transition. Right side. The trace anomaly in the confined 
phase measured on lattices with various lattice spacings and the continuum extrapolation (yellow band). The 
dashed line corresponds to the gluon gas contribution, estimated from the twelve lightest glueballs. 



Zooming in to the low temperure region we can explore the thermodynamics of the confined 
phase. To find out to what extent glueballs dominate, we calculated the trace anomaly contribution 



of the first twelve glueballs in Ref. [ |16| ] and plotted this together with our lattice results in the right 
side of figure. It has been suggested in Ref. [17] that the apparent deficit between the lattice data 
might be explained if we allow a temperature dependence for the glueballs. This dependence has 



already been determined for the O-i-i- and 2++ states [18]. This point was raised when only a couple 
of inprecise simulation points existed below from Panero's data set [Q]. Our data is in fairly 
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good agreement with this scenario, however, further studies are needed to understand the question 
in more detail. 

The strong non-perturbative effects in the equation of state can be emphasized by plotting 
(e — 3p)/r^ instead of normalizing it to (left side of figure For dimensional reasons, any 
finite order perturbative formula will only contribute logarithmic corrections to the p{T) ~ 
Stefan-Boltzmann law. Now plotting {e — 3p) /T"^ we expect to see a curve ~ times logarithmic 
corrections. What lattice data, actually, does show, is an approximately linear section up to ~ 
4Tc temperature, which then connects to the perturbative estimate. This linear segment was more 
striking in the Bielefeld lattice data |^, since the somewhat bigger errors hid a more complex 
structure. Moreover the latter data set ended at a temperature where the non-perturbative behaviour 
was still dominant. The observed non-perturbative pattern induced speculations on a "fuzzy" bag 



model [ |19| | and potential explanations in terms of a dimension-2 gluon condensate emerged [ ]20| , 
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Figure 2: Left side. Our results for the normalized trace anomaly multiplied by T'^/T^ for Nt = 5,6,7 and 
8 (red, green and blue dots, respectively). Also plotted are lattice results of [^], perturbation theory [Ql 
and the HTL approach |jl|]. Right side. Trace anomaly with various volumes at one of our lattice resolutions 
Nt = 5. Unless the box is very small (e.g. r — 4), there is no significant difference whether or not the box 
size allows contributions from the inverse scale. 

In this work we do not go into the viability of various explanations to the apparent ~ 
behaviour of {e — 3p), but we identify it as the dominant non-perturbative effect in the deconfined 
phase. Its effect reduces at high temperatures and becomes unnoticeable from the lattice equation 
of state, independently whether or not the lattice volume accommodates the inverse scale. One 
way of discussing the relevance of the Tc scale in the dynamics is to compare the trace anomaly at 
various volumes, as we do in the right side of figure ^for one lattice resolution. The "standard" 
aspect ratio r = 4 gives somewhat smaller values, but beyond r > 6 our simulation is not capable 
of resolving the difference. 

We summarize our findings as i) the large volume lattice trace anomaly data shows qualitative 
(and as we find using the fitted order coefficient, also quantitative, see later) agreement with 
the perturbative results from T > lOTc, and ii) we see no deviation between results from various 
volumes (with r > 6), moreover iii) the dominant non-perturbative contribution loses significance 
as ~ 1/r^. These considerations suggest that - even if the lattice volumes are ever shrinking as 
the temperature is increased - our results are able to describe the physical trace anomaly (and its 
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integral, the thermodynamic potential) within the error bars shown. Of course, this assumes that all 
perturbatively relevant scales are properly accounted for. 

Even in the small-volume simulations presented in this paper the lattices do accommodate the 
ultrasoft scale l/{g^T). This encourages us to use these very high temperature lattice data to extract 
some information on the ultrasoft physics. At high enough temperatures, where the &{g^) order is 
known to give a small correction to the iff {g^ ) order one can fit the g^ coefficient for the perturbative 
result, and the preferred renomalization scale /^htl for the HTL result. For the latter we find that 
our lattice data prefers /iHTL/(27rr) = 1.75(2)(6)(50), with the numbers in the parentheses from 
left to right representing the statistical error, the error coming from the lattice scale and that from 
the variation of the fit interval. Repeating the fit in Ref. we get for the sixth order coefficient 
qc = —3526(4) (55) (30), in the same notation for the errors. We used these data in the left and right 
side of Fig' |3|to compare our small volume simulations with the theoretical results. 




Figure 3: Left side. The continuum limit obtained from the lattice results (red band), compared to fitted 
perturbation theory. We fit the coefficient (gray dashed-dotted line), and the non-perturbative contribution 
in two different forms (black dashed line and green dotted line). Right side. The trace anomaly measured on 
two different spatial volumes (red and blue bands), compared to the NLO and NNLO HTL expansion with 
varied renormalization scale 0.5 < pL\iji^/2nT < 2 (green and gray shaded regions). The dashed-dotted line 
represents the expansion with the fitted scale. 

In the figure we also plot the perturbative formula plus a nonperturbative ~ contribution 
using two simple functional forms anp/r^ and (flnp + ^nplog(r/rc))/r^. Using the former ansatz 
we obtain a^p = 0.879(2)(40), and the latter anp = 1.371(1)(50), bnp = -0.618(2)(4). We note 
that these parameters are rather sensitive to the variation of the lower endpoint of the fit interval, but 
the relatively small ;^^/dof ~ 4 indicates (as it can be seen in the figure) that the low-temperature 
region is well described by the latter ansatz. Using thermodynamic relations the fitted perturbative 
formulae for the pressure and the energy density are also straightforward to calculate. In figure || 
we compare our small volume results to the so obtained predictions. Similar comparisons can be 
made for the case of the entropy density, too. 

Tabulated data and the details of the present work will follow in a separate, longer publication. 
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Figure 4: The normalized pressure (left side) and energy density (right side) as measured on our small 
volume boxes. A comparison is shown to various fitted perturbative functions, in the same notation as in the 
left side of the previous figure. 
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